Constraints on flavor-dependent long range forces 
from solar neutrinos and KamLAND 



Abhijit Bandyopadhyay^ * , Amol Dighe^^, and Anjan S. Joshipura^-^ 

1. Tata Institute of Fundamental Research, 
Homi Bhabha Road, Mumhai 4OO 005, India 
. 2. Physical Research Laboratory, Ahmedabad 380 009, India 

O ■ Abstract 

<N : 

. Flavor-dependent long range (LR) leptonic forces, like those mediated by the L^. — or Lg — L^- 

I gauge bosons, constitute a minimal extension of the standard model that preserves its renormal- 

' izability. We study the impact of such interactions on the solar neutrino oscillations when the 

. interaction range Rlr is much larger than the Earth-Sun distance. The LR potential can domi- 

> : 

, nate over the standard charged current potential inside the Sun in spite of strong constraints on 

CN ' the coupling a of the LR force coming from the atmospheric neutrino data and laboratory search 

o ; 

■ for new forces. We demonstrate that the solar and atmospheric neutrino mass scales do not get 

, trivially decoupled even if ^13 is vanishingly small. In addition, for a > 10"^'^ and normal hierar- 

CL|| chy, resonant enhancement of ^13 results in nontrivial energy dependent effects on the Ve survival 

probability. We perform a complete three generation analysis, and obtain constraints on a through 
• • . a global fit to the solar neutrino and KamLAND data. We get the 3a limits < 3.4 x 10~^^ 

^ I and Oer < 2.5 x 10^^^ when Rlr is much smaller than our distance from the galactic center. With 

' larger Rlr, the collective LR potential due to all the electrons in the galaxy becomes significant 

and the constraints on a become stronger by upto two orders of magnitude. 
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I. INTRODUCTION 



The standard electroweak model is now a well-established theory but it is believed to be 
incomplete and one expects some physics beyond the standard model (SM) to exist. Most 
extensions of SM postulate new physics at scales higher than the electroweak scales starting 
from TeV to the grand unification or Planck scale. There however exists an interesting 



iis_riiay arise 
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possibility that new physics may exist at scales below the electroweak scale. T 
from the existence of exactly or nearly massless gauge or Higgs bosons 
have remained invisible because of their very feeble couplings to the known matter. Various 
scenarios involving new physics at low energy and their possible signatures 
have been studied. 

Gauged extension of the SM is one possible scenario with new physics below the elec- 
troweak scale. Such a possibility is strongly constrained theoretically from the renor- 
malizability, however there exist Q three possible U{l)x gauge extensions of the stan- 
dard model which are anomaly free with minimal matter content. These correspond to 
X = Le — L^, Le — Lt, — Lt-. The extra gauge boson corresponding to U{l)x may 
not have been discovered if it is very heavy or if it is (nearly) massless but couples to the 
matter very weakly. The former possibility is analyzed in 1^1 . The latter possibil ity, first 
suggested in |8j , is strongly constrained by the search for the long range (LR) forces (ill . Il2 1 . 

Unlike the gravitational force, the U{l)x induced force couples only to the electron (and 
neutrino) density inside a massive object. As a consequence, the resulting acceleration 
experienced by an object depends on its leptonic content and mass. Such forces that violate 
the equivalence principle are strongly constrained. In case of the force with a range of ~AU, 
the most stringent bound comes from lunar ranging 
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121 which measures the differential 



acceleration of the Earth and moon towards the Sun. If a denotes the strength of the long 
range potential then these experiments inaply a < 3.4 x 10~^^ {2a) for a range A > lO^^cm. 

The flavor-dependent long range force [l] induced for example by Lg — L^^r (3, Isl can still 
influence neutrino oscillation in spite of such strong constraints on a. This happens because 
(i) the X-charge of the electron flavor is opposite to that of muon or tau flavor, so that 
these twox flavors propagate differently in matter and (ii) the large number of electrons (e.g. 
inside the Sun) and the long range of interaction compensates for the smallness of coupling 
and gives rise to a significant potential. For example, the electrons inside the Sun generate 
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a potential K/3 at the Earth surface given by 



V^%R^,) = ae,^ ^ 1.3 X 10-i^eV (^) , (1) 
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where ae/3 = ^ corresponds to the gauge couphng of Lg — Lp {[3 = yU, r) symmetry which 
we will sometimes collectively refer to as a. Here iV® ~ 10^'^ is the total number of electrons 
and Res is the Sun-Earth distance ~ 7.6 x 10^^ GeV ^. This is to be 
compared with the typical value of ^m? / E ~ 10^^^ eV for the atmospheric neutrinos. It 
follows that Vej3 can induce significant corrections to neutrino oscillations at the Earth even 
for a ~ 10-5°. 

One can define a parameter 

which measures the effect of the long range force in any given neutrino oscillation experiment. 
The bound on a from Q] implies that ^ < 750 in atmospheric or a typical long base 
line experiment, while < 35 for the typical parameters of the KamLAND experiment. 
Relatively large values of ^ tend to suppress the atmospheric neutrino oscillations. The 
observed oscillations can then be used to put a stronger constraint on a which were analyzed 
in ^. One finds the improved 90% C. L. bound 

a^f, < 5.5 X 10^52 ^ < Q_4 X iQ-^^ (3) 

in case of the Lg — L^^r symmetry respectively. 

With the improved bound on a given in Q, the value of ^ for KamLAND becomes 
rather small: ^ < 0.06. So one expects the KamLAND results to be infiuenced by the LR 
interactions to a very small extent. However, the potential Vep at the surface of the Sun is 
V^f^iro) ^ 2.8 X 10-^(ae/3/10-5°) eV, which may be compared with the MSW contribution 
Vcc ~ 6.0 X 10~^^ eV at r = 0.05 r©. Therefore one expects the long range potential to 



change or disturb the MSW LMA solution of the solar neutrino problem jlj] . Note that the 
effects on the solar and KamLAND experiments are qualitatively different, since KamLAND 
only probes the potential at the Earth given in eq. (0) while the solar neutrinos experience 
a long range potential that varies with the distance from the center of the Sun. It is thus 
important to do a combined analysis of these two experiments. 

The aim of this paper is to discuss new physical effects associated with this force and 
also make a quantitative analysis of the combined solar and KamLAND data to obtain 



a bound on a. It turns out that the long range potential produces physically interesting 
and quantitatively significant effects which can be used to constrain its strength. The bound 
obtained on a is more stringent than that obtained Q from the atmospheric results alone by 
more than an order of magnitude. If Rlr ^ -Rgai where i?gai is our distance from the galactic 
center ~ 10 kpc, the constraints become even stronger by upto two orders of magnitude. 

The plan of the paper is as follows. In Sec. |n] we present our basic formalism where 
we describe the main features of the LR potential inside and outside the Sun. In Sec. IIIH 
we present an analytic discussion of our results on neutrino masses, mixing angles and the 
resonances they undergo in case of the Le — symmetry. The corresponding analysis for 
the Le — is similar and the relevant analytic expressions are given in the appendix |XJ 
Sec. IIVI analyzes the KamLAND and the solar neutrino data numerically to obtain bounds 
on a for Rlr ^ -Rgai- The case Rlr^ Rgai is analyzed in Sec. El A summary of the results 
is given in Sec. IVII In addition, appendix |Bl gives a brief discussion of the impact of the LR 
potential on the neutrinos from a core collapse supernova. 



II. FORMALISM 



We consider the standard electroweak model with its minimal fermionic content but 
assume the presence of an additional gauged U{l)x symmetry. The cancellation of anomalies 
requires X = L^ — L^, Le — Lr or L^ — Lr [1]. The last symmetry does not play any significant 
role in the solar neutrino oscillations because of the absence of muons or tau leptons inside 
the Sun (or Earth). We will therefore concentrate on the first two and the couplings of the 
mediating vector bosons. The value of a is positive in this case. 

The observed neutrino oscillations imply that the U{l)x gauge symmetry cannot be an 
exact symmetry in nature. This is easy to argue. If it were exact, then the effective five 
dimensional neutrino mass operator following from any mechanism (e.g. seesaw) would be 
invariant under it. Consider the case of Lg — L^. Invariance under this dictates the following 
structure for the effective neutrino mass matrix: 
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e/1 












V 



(4) 



m, 

This structure implies a Dirac and a Majorana neutrino which remain unmixed and therefore 



cannot give any neutrino oscillations. Thus — needs to be broken. The symmetry 
breaking scale required to generate the solar scale Amf2 would be Am^g/'^^e^t- With Am^g ~ 
10~^ eV^ corresponding to the solar mass difference and ~ 0.1 eV corresponding to the 
degenerate neutrino mass, Am^g/^e/x is required to be at least 10~^ eV. A similar conclusion 
also holds in the case of the unbroken Lg — Lj. symmetry. 

The size of the U{l)x breaking as required above can be consistent with a nearly massless 
gauge boson since the corresponding coupling g in this case is required to be very small 
(< 10^^^) from the search of the long range forces 
the coupling also ensures that the relatively large U{l)x breaking in the neutrino sector is 
consistent with a very light gauge boson. In fact, a Higgs vacuum expectation value of a few 
GeV can lead to a gauge boson corresponding to the Earth-Sun range with g ~ 10~^^ — 10^^^ 
and can imply a relatively large neutrino splitting j8| . 

The most significant effect of the light gauge bosons would be in the solar neutrino 
oscillations. The coupling of the solar electrons to the L^. — L^_^ gauge bosons would generate 
a long range potential. If ne(r) denotes the spherically symmetric electron number density 
inside the Sun then the long range potential is given by 

Ve%r < ro) = Ana^p / ^ / r"'n4r")dr" . (5) 

Jr ^ Jo 

Outside the Sun, 

V;0(r > Tq) = r r"W{r")dr" = ^Nf . (6) 

r Jq r 

The approximate profile 

rieir) ~ 245 Na 10-^°-^^ ^^Z'^®) cm'^^ (7) 

of the solar density implies that is a monotonically decreasing function, which is 
inversely proportional to r when outside the Sun. This behavior is shown in Fig. ^ which is 
obtained using the actual electron density profile in the Sun. It is seen that dominates 
over the MSW potential Vcc inside the Sun for a > 10~^^. Moreover, it does not abruptly 
go to zero outside the Sun like Vcc, but decreases inversely with r, ultimately reaching the 
value given in eq. ((H) at the surface of the Earth. When Ve/3 > Vcc inside the Sun, the 
resonance is shifted outwards (sometimes even outside the Sun) and its adiabaticity may be 
affected. 




The required smallness of 
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FIG. 1: Comparison of the MSW potential Vcc and the LR potential due to the solar electrons 
from the solar core all the way to the Earth (r/r© ~ 215) and beyond. The {Am? /2E) values 
corresponding to = 10 MeV are also shown. 

The contribution of the electrons inside the Earth can be calculated in a similar 
fashion. Roughly, one finds that at the Earth surface, 

where Mq{Me) and Res{RE) respectively refer to the mass of the Sun (Earth) and the 
Earth-Sun distance (radius of the Earth). Thus the solar long range potential dominates 
over the terrestrial contribution and we will neglect the latter. As long a.s Rlr ^ Rgai, this 
is the dominant potential affecting the propagation of solar neutrinos. 

When Rlji > Rgai, the collective potential due to all the electrons in the galaxy may 
become significant. The mass of the Milky Way is (0.6 - 3.0) xlO^^ solar masses, which is 
mostly concentrated in the center of the galaxy. The baryonic contribution to the galactic 
mass may be estimated to be 0(10%). The center of the galaxy is ~ 10 kpc away from the 
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< = ^ «e/3 ^ , (9) 



Sun. We denote the galactic contribution to the potential Ve^ as 



gal 

where A^°gai is taken to be 10^^ A*"® and R^^^^ to be 10 kpc. The net LR potential is Vep = 
V^ + Vefi^- The parameter b takes care of our ignorance about the distribution of the baryonic 
mass in our galaxy. With R^fi > -Rgai, we expect 0.05 < b < 1. The value of b may be 
smaller if Rui is smaller. Clearly, 6 = would have the same effect as <C -Rgai- With 
6 7^ 0, the constraints on a become stronger, as will be demonstrated in Sec. |^ 

The screening length due to the antineutrinos present in the cosmic neutrino background 
is a few hundred kpc for rriu ~ 0.1 eV Ql- Therefore, for the galactic scale, the screening 
plays no significant role. Over the Sun-Earth distance, even the possible local screening 



effects would be too small to have any effect 

In addition to the altered resonance structure inside and outside the Sun, the mixing 
angles at the Earth also differ from the corresponding vacuum values, with the result that 
both the solar and the KamLAND neutrinos get affected by the LR potential. An important 
point to note is that this potential gives unequal contributions to two flavors (e and ^ or r) 
simultaneously unlike in case of the charged current which contributes only to the electrons. 
The third flavor gets no contribution. As a consequence of this, the inclusion of three 
generations in the solar analysis becomes necessary. 

The appropriate Hamiltonian in the flavor basis describing the neutrino propagation can 
be written as 

Hf = -R23(^23)-Rl3(^13)-Rl2(^12)-f^0-Rl2(^12)-Rl3(^13)-R23(^23) + V , (lO) 

where Hq refers to the effective Hamiltonian in the mass basis, and Ri/s are the rotation ma- 



trices in the i-j plane. Since the absolute masses of neutrinos play no part in the oscillation 
phenomena, we can tak( 
respectively, leading to 



phenomena, we can take the neutrino mass eigenvalues in vacuum to be 0, a/A?t2|^, a/A??t|^ 



ifo = Diag(0,A2i,A32) , (11) 

where A21 = A'm2i/{2E) and A32 = Am32/(2£'). The rotation angles 623 and 612 are 
the vacuum mixing angles describing the atmospheric and the solar neutrino oscillations 
respectively, whereas ^13 is the third "Chooz" mixing angle. We have assumed that no CP 
violation enters into picture here. 



The matrix V in pOj) describes the combined contribution of the charge weak currents 
as well as the long range forces. Explicitly, 

V = Diag(F,e + Km, -Km, 0) . (12) 

The neutrino propagation is described by eq. (fTUI) . The corresponding antineutrino propa- 
gation is obtained by the replacement V —V. 



III. MASSES, MIXINGS AND RESONANCES OF SOLAR NEUTRINOS 



In order to analyze the propagation of solar neutrinos, we rewrite eq. (fTUI) explicitly as 



Hf = A 



32 



/ 2 \ 

XSi2+ He + yeti 3;Ci2Si2C23 + S13S23 — iCCi2Si2S23 — S13C23 

XC12S12C23 + S13S23 S23 + XC12C23 - 1/eM C23S23(1 - a;Ci2) 

\^ -a;Ci2Si2S23 - S13C23 C23S23(1 " Xcf^) C^g + XCI2SI3 J 



(13) 



where 



A 



X = 



21 



A 



0.03 , y,= 



Vrr 



32 



A, 



32 



2^Kc 

A"^i2 



Km _ 2^Km 



and Sij = sinOij, = cosOij. Since 613 is small (6^13 < 0.2 [l^), we have kept terms to only 
linear order in S13. 

Eq.()13|) can be diagonalized through the unitary matrix 



^32 



Am 



2 ' 

32 



(14) 



Um = -R23(^23m)-Rl3(^13m)-Rl2(^12m) 



(15) 



such that 



2E 



(16) 



The smallness of x and S13 can be used to approximately determine the matter dependent 
mixing angles of Um to the leading order in these parameters. The angle 6'23m follows from 
the lower right 2x2 block in (|T3|l : 



tan 26'' 



sin 26*23(1 — xc 



23m 



12) 



COS 26*23(1 - XCI2) +yef. 



The subsequent diagonalization leads to 



tan 26* 



13m 



2 (a;si2Ci25' + S13C) 



C2 + xicl^S^ - SI2) -Vc- Ve^il + sin2 6; 



23m 



(17) 



(18) 



where S = sin(6'23m - 6*23) and C = cos(6'23m - 6^23) • 



As long as the denominator in eq. ()18p does not vanish (which happens only in a very 
narrow range of near i/efj, ~ 2/3), we can take 9i3m ~ 0{x,si3). Neglecting terms that 
are quadratic or higher order in 0{x, S13), the effective Hamiltonian in the new basis (after 
the 2-3 and 1-3 rotation) becomes 



^ Xsj^ + yc + Vepi XC12S12C - Si'iS ^ 

Hf. ^ A32 XC12S12C - + XCI^C^ - Vet^cj^m 

y + xcl^S^ - ye^,sl^ra J 
so that a 1-2 rotation through an angle 612m, given by 



(19) 



2e,2m - ,2^ , 2 J 'T"^ " n ^ 2. ^ ' (20) 

+ x{ci2C^ - 3(2) - Vc - ye^,{l + COS^ 023m) 

completes the diagonalization. The neutrino masses are given as 

A32^ [x{cl2C^ + S^) +y, + sin^ ^^sm + 52 - D'/^] , 

A32^ [x{cj2C^ + S^) +y, + Sin^ ^^Sm + S^ + D'/^] , 



^Im 



m 



3m 



2A32E{C' + xci2S' - sin' 623^) , (21) 



where 

D=[S' + X{C12C' - S\2) - Z/e - ye,{l + COS^ ^23™)] ' + 4 (XS12C12C - 8^38? • (22) 

The above analytical results can be verified by the exact numerical results in Fig. |21where 
we show the angles and mf values in matter for different values of a for normal as well as 
inverted hierarchy. 

In this and the next section, we analyze the case Rlr <C -Rgah so that the potential 
is as shown in Fig. CI As is apparent from the figure, the maximum value of is given by 



a VIO MeV^ 

for the best fit values of the atmospheric parameters. Thus, at = 10 MeV, we have 
Ven ~ 0.1 for a = 10~^^. The left column in Fig. |21then corresponds to the range of a where 
{yeii)max < 1, and the right column corresponds to {ye^)max > 1- 

The propagation of solar neutrinos is qualitatively and quantitatively different depending 
on whether {ye^i)max is large enough to cause resonant enhancement of 9i3m in (fTH|) . The 
resonance occurs when a > 10~^^. We therefore consider the two cases a < 10^^^ and 
a > 10~^^ separately in the next two subsections. 
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FIG. 2: The angles and m? values in matter for solar neutrinos for E = 10 MeV, in the case 
Rlr ^ Rga.1- The mf values are correct up to an additive constant, so that only their relative 
values have a physical significance. 



A. For a < 10 



-52 



For a <^ 10~^^, we have ygfj. -C 1 and the atmospheric mixing angle gets only a small cor- 
rection from the matter effects. Writing ^23m = ^'23 + 56^23, we see that 6623 ~ —?/e/x sin 26^23/2 
and 6'23m remains close to its vacuum value, 623m ~ 7r/4. With a higher a value, the deviation 
6023 becomes appreciable and results in the reduction of sin^ '2923m as shown in Fig. |21 

The angle Oism has contributions from two sources: from finite 613 in vacuum as well as 
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the additional contribution from the term a;si2Ci25'. The latter is doubly suppressed because 
of the smallness of x as well as 5 ~ 862^, and can be neglected as long as S13 > xS . One 
may then take 6*13^ ~ C'lsis) since the resonant enhancement of 6*13^ anyway does not occur 
for a < 10"^^. 

In the limit 9i3m 0, the third mass eigenstate decouples and the scenario reduces to 2u 
mixing, as can be seen from eq. ()19|) . However, note that the effective matter potential is 

V12 ~ Vcc + V;^(l + cos' ^23m) , (24) 

and not V^c + 2Ve^ as would have been taken in a naive 2-generation analysis. Thus, the effect 
of the third neutrino and its mixing is inescapable here. However, it only appears through 
the factor (1 + c^^^) in eq. (j24p . and the mass of the third neutrino or the mass hierarchy is 
immaterial for the effective 2z/ analysis. This may be verified from the left column of fig. |21 
The most important effect of the LR potential is for the solar angle. Eq. ()20|) gives the 
resonance condition 

Am^i cos 2^12 [Kc + KM(l + cos2^23m)] , (25) 

which differs from the MSW condition by an addition of the term involving Vg^. For a > 
10"^'^, the Ve^ contribution dominates over Vcc and changes the MSW resonance picture 
significantly. The resonance is shifted away from the center as a increases. Eventually 
for some value of a the resonance gets shifted outside the Sun where its behavior is solely 
determined by the LR potential. For a > 10~^^, neutrinos with E = 10 MeV encounter 
resonance outside the Sun in case of the best fit values of the neutrino mass parameters 
obtained in the standard analysis jl^j ]. 

Addition of the Ve^(l + cos' 623m) term to Vcc makes the variation of the total potential 
smoother than the normal MSW potential with the result that the transition becomes more 
adiabatic than the corresponding case without the LR. In particular, when the resonance 
occurs outside the Sun then V^® ~ and the adiabaticity parameter at the resonance is 
given by 

_ Am^g sin' 26*12 1 dVu 



2E cos 2^12 



V12 dr 



^ aet^N® tan' 2^12(1 + cos' ^23) ~ 1.4 x 10^^ Oc^ . (26) 

The value of 7^ is independent of the neutrino masses, energy and position of the resonance 
and is solely determined by a and the vacuum mixing angles. For the standard values of 
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the latter, the resonance is found to be highly adiabatic: 7^ 1 for a > 10~^ . In general, 
if Pl{E) = exp{—TT'~fL/2) is the probability that uim and z/2m convert to each other while 
passing through the resonance, the net survival probability of is 

Pee{E) = (1 - Pl) COS^ O^p COS^ 9i2P COS^ 9i3E cos^ O^e 

+ Pl cos^ 6i3P sin^ Oup cos^ 6ise cos^ Que 

+ (1 - Pl) cos^ 6i3P sin^ 6'i2p cos^ 6'i3£; sin^ Oue 

+ Pl cos^ 6i3P cos^ 6'i2p cos^ 6'i3£; sin^ Que 

+ sin^ ^i3p sin^ OisE ■ (27) 

Here 6ijp and are the values of 6ijm at the neutrino production point and at the Earth 
respectively. The energy dependence of Pl as well as all the angles is implicit. Note that 
since 0i3p,9ise ~ 0{9is), the last term may be neglected if we neglect terms of 0{9f^) or 
smaller. 



B. For a>10"^^ 

For a ~ 10~^^, the value of is large enough so that sin^ 26'23m gets unacceptably 
suppressed through eq. (fTTj). This also suppresses the atmospheric neutrino flux and results 
in the bounds on a discussed in j^. 

For solar neutrinos, the i'im-^2m resonance as described in the previous section occurs, 
but in addition the angle 6*13^ gets resonantly enhanced when 

+ xiS'cl^ - sl^) -y,- y^^il + sin^ ^23™) ~ . (28) 

This happens when ?/e/x ~ 2/3. The sign of also needs to be positive, so the resonance 
occurs only for normal hierarchy. For the inverted hierarchy, there is no resonance for 
and eq. (f?fj) gives the correct expression for their survival probability. 

In the resonance region, the effective Hamiltonian matrix (fTn|) can no longer be diago- 
nalized through the simple procedure described in the beginning of Sec. IIIH and the mixing 
angles have to be computed numerically. However, this happens only in a small range of ye^ 
around ~ 2/3: the width of the resonance region may be estimated to be 6yefj. ~ (2/3)si3. 
The expressions ()20| ) - (j^^ are valid everywhere outside this region. 
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FIG. 3: dependence of Ph for various values of a ioi E = 10 MeV in the case Rlr *C i?gai- 

The 6*13^ enhancement corresponds to the i'2m-^3m level crossing, with an effective po- 
tential 

V23 = + + sin^ d23m) . (29) 

When the hierarchy is normal, only a fraction of the that are produced mainly as i>2m 
inside the Sun survive the z/2m-t'3m resonance. The adiabaticity at this resonance, which 
strongly depends on 6*13, affects the net survival probability of z/g: 

Pee{E) = COS^ 6i3P COS^ 6i2P COS^ 6i3E COS^ di2E 

+ (1 - Ph) cos^ 6i3p sin^ 612P cos^ 613E sin^ O^e 

+ Ph sin^ 6i3p cos^ 9i3e sin^ O^e 

+ {l-PH)sm^ei3psm^ei3E 

+ Ph cos^ 6i3p sin^ 9i2p sin^ 9i3e , (30) 

where Ph{E) is the probability that V2m converts to after traversing through this reso- 
nance. Here we have used the earlier result that for a > 10"^^ the vi^n ~ ^2m resonance is 
outside the Sun and is always adiabatic [see eq. (|^]. The energy dependence of Ph as well 
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as all the angles is implicit. 
The value of Ph is given by 



H 



exp 



2 



rrin 



(31) 



2E dei3m/dr ^, 

Clearly, if Ph ~ 0, the expression (jHUj) reduces to ^I7\ . and the results of the 2z/ analysis 
stay valid. In general Ph ~ at high values of 6^13. In Fig. El we show the 6*13 dependence 
of Ph for various values of a for E = 10 MeV. At a = 10~^^, the value of Ph > 0.1 for 
^13 < 0.08°, which is when the survival probability is affected significantly. For larger a, the 
value of Ph becomes significant for lower ^13 values. In the range where 0.1 < Ph < 0.9 
(the semi-adiabatic range), Ph is also highly energy dependent, as can be seen from (jHTj) . 

The analytic discussion above reveals that the LR potential makes important contribu- 
tion to the solar neutrino problem and a detailed numerical analysis is needed to obtain 
constraints on this potential. We turn to this analysis in the next section. 



IV. CONSTRAINTS FROM SOLAR NEUTRINOS AND KAMLAND 

To find the best fit values of the oscillation parameters and a from a statistical analysis of 
the experimental data, we employ the minimization technique with covariance approach 
for the errors. For analysis of the total event rate data from all the experiments, the 
function is defined as 

iVexpt 

where Pf (^ = th or expt) denotes the total event rate for the z*^ experiment. Both the 
theoretical and experimental values of the fitted quantities are normalized relative to the 
standard solar model (SSM) predictions. The error matrix {(7if°'^^^Y contains the experi- 
mental and theoretical uncertainties along with their correlations. Theoretical uncertainties 
include the uncertainties in the capture cross sections, which are uncorrelated between dif- 
ferent experiments and the astrophysical uncertainties from the SSM predictions which are 
correlated between different experiments. The correlations are being evaluated using the 
procedure of IgI ]. 

For the analysis of any spectral data (recoil energy spectra or zenith angle spectra), the 
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is defined as 

A'bins 

x^pec = E - ^r^') ['^^^rr] ' (^f - ' (33) 

where Sf (^=th or expt) is the number of events in the i^^ bin of the spectrum. The 
error matrix for the spectral data includes the statistical error, correlated and 

uncorrelated systematic errors in the different bins and the error due to the calculation of 
the neutrino energy spectrum from SSM. 

For a global analysis of the solar data - rates from CI, Ga experiment, spectral data 
from SuperKamiokande (SK) and Sudbury Neutrino Observatory (SNO): both D2O and 
salt phase, and KamLAND data - the relevant is given by 

X =X CljGa rates + x SK spec + x SNO spec + x KamLAND • 

(34) 

Note that only solar neutrino observations would not have been able to put strong con- 
straints on a: as long as there is no i'2m-'^3m level crossing and the i'im-'^2m resonance is 
adiabatic, the net z/g survival probability (j27j) is a function of Am^g/*^ ^ large a, i.e. for 
Vcc << Vefj,. As a result, one can fit the solar data by increasing the value of Am^2 when a 
is increased and a strong bound on a would not follow. 

However, the data from KamLAND restricts Amg^ to a very small range and plays a 
crucial role in constraining a. We use eqs. (I27|l and pO|l for the survival probability of solar 
neutrinos. The z7g survival probability in KamLAND is given by 



ee 
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pis^ = 1 - cos^ e^^ 

,2 



■ 2oQ • 2 / ^"^21-^ 

sm 26^12 sm ' 



r „ / Am: 

+ sin^ 26*13 sin^ 6*12 



sm^2^i3sm^(^ 



(35) 



\ AE J \ AE 

where all the mass squared differences and the angles are measured at the Earth for antineu- 
trinos. Note that for antineutrinos, the sign of Vcc as well as Vg/i is reversed with respect to 
the neutrinos. 

In Fig. U we show the Ax^ values as a function of the parameter a for various 6*13 values. 
The best fit values for the solar parameters are always observed to lie in the LMA range with 
vanishing giving the best fit. For a < 10~^^, the value of is minimum for ^^13 = 0°, 
which is consistent with the observation that 6*13 = 0° also gives the best fit to the solar and 
KamLAND data when the LR forces are not taken into account When a > 10^^^, a 
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FIG. 4: Ax^ = x^ic^e^) — X^i^^e^i = 0) values for different ^13 values, in the case Rlr <C ^gal- 
strong energy dependence in the survival probability is introduced for 6 < 0.8° through Pfj, 
so that the values for extremely low ^13 values become large. In this region, the lowest 

is found to be at values of ^13 that are small, but still keep Ph ~ 0. We have shown 
corresponding to such a 613 in the figure. 

The bounds on a should therefore be, strictly speaking, 6'13-dependent. However, the 
region a > 10~^^, where the 6*13 dependence from Ph starts coming into picture, is excluded 
to more than 3cr as can be seen from Fig. EJ Therefore the constraints on a by using 6*13 = 0° 
are the most conservative ones, and we quote the upper bounds on a obtained by taking 
^13 = 0°. These limits are shown in Fig. El the 3a limit corresponding to the one-parameter 
fit is 

aef, < 3.4 X 10"^^ . (36) 
The corresponding limit in the — Lj. case (see Appendix ^ is 

Oer < 2.5 X 10"^^ . (37) 

The bounds are independent of whether the neutrino mass hierarchy is normal or inverted. 
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FIG. 5: Ax^ values and limits for the Le — as well as — symmetry, with ^13 = 0°. The 
case Rlji <^ -Rgai is represented by 6 = and higher b values correspond to larger contributions 
from galactic electrons (see Sec. IVjl . 

V. THE LONG RANGE POTENTIAL DUE TO GALACTIC ELECTRONS 

The collective contribution of all the electrons in the galaxy to the LR potential in 
the solar system may be parametrized in general as given in eq. 0. The net potential 
Ve^ = V^® + V^^^ is shown in Fig. IHlfor various values of b and a. Clearly, larger the value of 
6 or a, larger the value of Vg^. Also, note that the value of Vg^ near the earth is approximately 
the same as V^^^, since V^® keeps on decreasing as one travels towards the Earth, whereas 
V^^^ is a constant over the scale of the solar system. 

With our understanding of the effects of the LR potential on neutrino masses, mixings 
and resonances obtained in Sec. IIIH the following observations may be made: 

(i) For V^^^ ^ Am^j/{2E), there is no MSW resonance that is essential for a good fit to 
the solar neutrino data. Therefore, larger values of h and a {ha > 10"^'^) are expected to be 
ruled out from the global fit. 

(ii) For V^^^ <C /S.m1^^^J {2E) , the matter potential Vcc dominates inside the Sun, and the 
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FIG. 6: The net potential Vef, = V® + if^f' for various values of b and a. 

standard picture of the neutrino flavor conversions inside the Sun is not affected. Therefore, 
smaller values of h and a {ha < 10~^^) should be allowed. 

(iii) For the intermediate values of Vg^^\ the situation depends strongly on whether the 
potential proflle near the MSW resonance is dominated by Vcc or K^T^- ■'■^ former case, 
the resonance is adiabatic for > 5 GeV and only partially adiabatic for lower energies, 
which gives a good flt to the data. In the latter case, however, the resonance tends to be 
adiabatic even for low energies, so that the radiochemical data will disfavor the solution. 

The value of h is expected to be in the range 0.05 < 6 < 0.5 (see Sec. 10]) • The Ax^ values 
as a function of a for 6 = {i.e. Rlr -Rgai), b = 0.1 and 6=1 are shown in Fig. The 
3cr constraints for — are 

Oe/. < 2.9 X 10"^^ (6 = 0.1) , < 2.6 X 10-^^ (& = 1) (38) 

and for L^. — L^, they are 

fter < 2.3 X 10"^^ {b = 0.1) , aer < 2.1 X 10"^^ (&=!)• (39) 

Clearly, the constraints get stronger as b increases. The most conservative constraints are 
therefore with 6 = 0, as calculated in Sec. IIVI 
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VI. SUMMARY AND CONCLUSIONS 



Flavor-dependent long range leptonic forces, like those mediated by the Le — L^ or Le — L-r 
gauge bosons, constitute a minimal extension of the standard model that preserves its renor- 
malizability. The flavor dependent potentials produced by these forces influence neutrino 
oscillations. The effects of these are quite significant in spite of the very strong constraints 
on the couplings of such forces from astronomical observations or Eotvos type laboratory 
experiments. We have performed a detailed study of specific effects of these forces in the 
solar neutrino and the KamLAND experiments. 

It was found that the new forces change the standard MSW picture in a quahtatively 
different way which ultimately results in a strong bound on the couplings of these forces. 
We have developed a detailed formahsm to describe these effects and have used it to obtain 
bounds on the couplings from the statistical analysis of the experimental data. It was shown 
that the mixing among all three generations needs to be taken into account because of the 
fact that the Lg — L^^t gauge bosons couple to two out of three flavors at a time. The changes 
which result in the MSW analysis were studied both analytically as well as numerically in 
the case Rlr ^ Rgah when the galactic electron contribution to the LR potential may be 
neglected compared with the solar electron contribution. 

A qualitatively new effect studied in detail is the possible resonant enhancement of ^^13. In 
the standard MSW picture, a non-zero but small 6*13 can only give sub-leading corrections. 
In contrast, the long range potential can resonantly amplify ^13 if a ^ 10~^^ and the neutrino 
mass hierarchy is normal. The global analysis of the solar data however constrains a < 10~^^. 
As a result, the resonance enhancement of ^13 does not take place in the solar case. But this 
resonance effect can play an important role in other environments, e.g. inside a supernova. 
See Appendix B for details. 

A global analysis of all the solar neutrino and KamLAND data was performed to 
constrain the coupling a. The solar data alone are found to be inadequate in constraining 
a: one could always fit these data by appropriate change in Ami2 compared to the standard 
LMA values. This does not remain true when the KamLAND results are included. A 
significant bound on a is obtained by combining the solar and KamLAND results. The 
conservative 3cr bounds follow when 613 = 0: 

CKe^ < 3.4 X 10-^^ , aer < 2.5 X 10-^^ . (40) 
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These bounds are stronger by more than one order of magnitude than the ones in eq. Q 
following from the analysis of the atmospheric neutrino data. 

A much stronger bound on a, namely a < 6.4 x 10~^^, was quoted in purely from the 
solar neutrino results. This was not based on the detailed statistical analysis as presented 
here, but was obtained under the assumption that even in the presence of a the Am^2 
9i2 should lie within their 95% C. L. range obtained in the standard LMA solution. This 
assumption need not a priori be true. In fact as discussed here, the solar neutrino results by 
themselves cannot be used to constrain a, so the detailed analysis as done here is required. 
An analysis similar to ours has recently been carried out in [l2| which has reported bounds 
on couplings of the vector and non-vector long range forces. The resulting bound on a 
in the former case is similar to ours. However, they assume one mass scale dominance, 
neglecting the mass eigenstate z/3 altogether. As shown in this paper, in the presence of the 
LR potential, 1/3 affects the solar neutrino survival probability significantly even when ^13 is 
vanishingly small. Moreover, the galactic contribution has not been included in the analysis 
of ^1 even when the range of the force is more than our distance from the galactic center. 

When Rlh > -Rgai, the collective contribution of all the electrons in the galaxy to the 
LR potential becomes significant. This gives more stringent constraints on the value of a, 
which also depend on the distribution of baryonic mass within the galaxy. We parametrize 
our ignorance about this with a parameter h (expected to lie between 0.05 and 1 with 
conservative estimates) and perform global fits to constrain a for fixed h values. We obtain 
tte^ < 2.9 X 10"^'* for h = 0.1 and Oe/x < 2.6 x 10~^^ for 6 = 1 in the Le — case. In the 
Le — Lr case, one gets < 2.3 x 10"^^ for b = 0.1 and < 2.1 x 10^^^ for b = 1. Clearly, 
the constraints become stronger as the galactic electron contribution, or the range of the 
potential, increases. 

The strength of the LR forces increases with the electronic content of the source and there- 
fore their effects are expected to be much stronger for supernova neutrinos. As discussed in 
Appendix B, the conventional flavor conversions of the supernova neutrinos changes signif- 
icantly in this case even for a ~ 10^^^. In particular, the LR induced resonance remains 
adiabatic for very low values of ^13 and the Earth matter effects may be absent. Also, the 
shock wave effects on the neutrino spectra may be absent for t < 10 s, which is when the 
neutrino flux is significant. On the other hand, the observation of any of these effects may 
be used to improve the bound on a at the level of 10~^^, even when the galactic contribution 
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to the LR forces is small. 

While the existence of LR forces may be regarded as a theoretically allowed specula- 
tion at this stage, it is quite remarkable that these forces, if they exist, strongly influence 
the atmospheric and solar neutrino oscillations. They would also effect the long baseline 
experiments which can provide additional constraints on a. 

We have concentrated on bounds on the gauge coupling a of the LR forces. In principle, 
the gauge symmetry allows mixing between the Lg — L^^,- gauge boson X and the ordinary 
hypercharge gauge boson B in their kinetic energy terms. This mixing would lead to mixing 
between the X boson and the photon and would lead to a flavor dependent infinite range 
potential even if the X boson has a finite mass. The strength of this force will be governed 
by an independent mixing parameter Q times the electromagnetic coupling a. Based on the 
present analysis, we expect this quantity to obey the same constraint as obeyed by OL^^^r in 
case of the infinite range potential. 
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APPENDIX A: CONSTRAINTS ON THE Le - GAUGE BOSON COUPLING 

The analysis of Lg — gauge bosons can be carried out in an analogous manner. The 
potential in the fiavor basis becomes 

y = Diag(Ke + V;„0,-Kr) (Al) 
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and the relevant expressions for the mixing angles in matter (to the leading order in x and 
S13) are: 



sin 26*93(1 — xcir, I , . , 

tan2^23m ~ 2 , ' A2 

cos26'23(l - xcf^) - y, 



er 



tan 2^13 ^ '2{xsi2Ci2S + gisC) ^^^^ 

C2 + a;(ci25'2 - S12) - 2/c - 2/6^(1 + ^os^ 6*23™) 

- ,2^ < . r^^'fT^ ~ ^ - .a ^ ■ 

+ x(cf2C^ - Sfa) - 1/c - 1/6^(1 + t'23m) 

The resonance structure is similar to that in the Lg — case. The limits on the coupling of 
such gauge bosons are shown in Fig. |3| The bounds are independent of whether the neutrino 
mass hierarchy is normal or inverted, like in the case of L^. — L^. 



APPENDIX B: EFFECT ON A CORE COLLAPSE SUPERNOVA 

The bounds on the LR forces that we obtained from the atmospheric, solar and Kam- 
LAND experiments are of the order a ~ 10^^^, when the galactic contribution to the LR 
forces is small. Although these bounds seem very stringent, even such a small strength of 
LR forces can potentially give rise to significant effects in the neutrino spectra from a core 
collapse supernova. The spectra of Ve and z/g from the SN have encoded information about 
the primary neutrino fluxes and neutrino mixing parameters and they can even show 
signatures of the passage of the shock wave through the mantle Note that all the 

above analyses have been carried out assuming that the collective flavor conversion effects 
caused by the neutrino-neutrino interactions are negligible compared to the conventional 
non-neutrino matter effects on neutrino propagation. If the collective effects happen to be 
strong, as claimed in [2^, our estimations in this section, as well as most of the SN flavor 
conversion analyses till now need to be reexamined. 

In Fig. we show a typical profile j2l| of the MSW potential Vcc inside a SN as well as the 
profile of the LR potential Ve^ for two values of a that are allowed with the constraints found 
in the conservative scenario Rlr ^ -Rgai- Note that even with a as low as 10~^^, the LR 
potential Ve^ exceeds Vcc inside the star, and hence affects the dynamics of neutrino flavor 
conversions. The effects, which may be significant in the allowed range 10~^^ < a < 3 x 10~^^, 
will be as follows: 



(i) The positions of the H and L resonances 22], corresponding to Amg and Am^^^^ re- 
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FIG. 7: The MSW potential and the LR potential for a typical SN for two different a values that 
are allowed with the bounds obtained in this paper for Rlr ^ -Rgai- The SN model taken 21| is a 
star with a mass of ISMq and primordial metallicity equal to that of the Sun. 

spectively, are shifted away from the center of the star by a factor of up to one order of 
magnitude. 

(ii) If Vep, dominates over Vcc in the resonance region, the resonance is highly adiabatic, since 
the LR potential is in general smoother than the MSW potential. Therefore for larger a. 
values, both the H as well as L resonances are adiabatic for practically all values of Qxj,. The 
SN neutrino spectra then lose the ability to reveal any information about in the absence 
of any shock wave effects. For example, no Earth matter effects may be observed. 

(iii) The shock fronts will reach the resonances at late times, t > 10 s, when the neutrino 
flux has reduced a lot. As a result, the shock wave effects would be much harder to observe. 

(iv) On the other hand, if any effects of non adiabaticity, e.g. Earth matter effects or shock 
wave effects, are identified in the neutrino spectra, the bound on a can be improved by 
almost an order of magnitude, to a < 10^^^. Supernova neutrinos thus form the most 
sensitive probe for the LR forces, at least when their range is smaller than -Rgai. 
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(v) For Rlr > -Rgai, the constraints obtained from the SN observation are expected to 
be comparable to those found from the solar neutrinos and KamLAND, since it is the 
approximate condition (6 a Mgai/-Rgai) -C /\m1^-^^^ / {2E) that determines the allowed range 
of a, like in Sec. El 
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